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The current noise in a classical one-dimensional charge density wave system is studied in the weak 
pinning regime by solving the overdamped equation of motion numerically. At low temperatures 
and just above the zero temperature depinning threshold, the power spectrum of the current noise 
S(f) was found to scale with frequency / as S(f) ~ /~ 7 , where 7 ~ 1, suggesting the existence of 
flicker noise. Our result is in agreement with experimental findings for quasi-one-dimensional charge 
density wave systems and provides the first evidence of 1// behavior obtained from first principles. 

PACS numbers: 72.15.Nj 

INTRODUCTION 

The noise in any system is characterized in terms of the shape of its spectral power density S(f), which may be 
measured directly in experiments. If S(f) ~ /~ 7 and 7 w 1, then such a noise is referred to as 1// or flicker noise. 
Its existence has attracted the attention of researchers from various branches of the natural sciences for many years. 
As for physical systems, 1// noise has been observed as fluctuations in the currents of diodes, vacuum tubes and 
transistors, the resistance of carbon microphones, metallic thin films and semiconductors [HQ-El- the magnetization 
in spin glasses [j| and in many other systems. 

Despite a substantial amount of experimental and theoretical effort, a unique underlyingmechanism for the flicker 
noise remains unknown and still poses an open question as to its universal origins However, there is now 

increasing evidence that the non-trivial exponent, at least in electronic systems, may come from equilibrium dynamics 
in the presence of disorder || when energy barriers exceed the typical thermal energy. 

I// 7 noise in charge density wave (CDW) systems 01 nas been studied mainly in NbSe3 and TaS3 materials 0. 
The first experiment was done on a bulk NbSe3 sample by Richard et al. [ll| who found 7 ~ 0.8. Studying transport 
properties of the quasi-one-dimensional CDW material TaS3 at low temperatures, Zaitsev-Zotov observed that 
slightly above the depinning threshold of the driving electric field, the exponent 7 for the current noise is equal to 
7 « 1.2. 

It should be noted that a phenomenological model based on fluctuations in impurity pinning force due to deforma- 
tions of the sliding condensate was proposed to explain the broad band 1// noise in CDW systems 0|. However, 
neither theoretical nor numerical estimation of 7 has been provided so far. The aim of the present paper is to compute 
7 from first principles with the help of a one-dimensional classical model for CDWs |l4j . The current was obtained 
through numerical simulation of the overdamped equation of motion. The 1 // scaling is evaluated using the so-called 
Wavelet Transform Modulus Maxima (WTMM) method |15|. The exponent 7 was found to depend on T. At low 
temperatures (T < 0.1), in agreement with the experiments |l2(, we obtain 7 m 1.2 in the crossover regime. Exponent 
7 drops with increasing T and the "exact" l//-noisc is observed at T ss 0.3 where 7 becomes 1. This interesting 
result is indicative of the possible occurrence of 1// noise. At high temperatures 7 takes on the white noise value 0. 

Notably, the observed 7 « 1 is not related to the second order depinning transition behavior at T = 0. Due to the 
asymptotic uniqueness of the sliding state |lfij . this critical point dynamics scenario leads to the 'trivial' exponent 
7 k 2 [T7L llij. Additionally, the observed 'flicker' noise behavior 7 k 1 gains on its scaling range with increased 
distance to the critical point of the second order depin ning transition. 

Based on unusual current -voltage characteristics [l^, Oil |2(J, Um j Zaitsev-Zotov suggested that at low temperatures 
the quantum creep dynamics may play an important role and proposed the crossover from classical to quantum creep 
regime as an alternative explanation for experimental results. The strength of the quantum fluctuations in ID CDW 
systems can be estimated by a dimensionless parameter K which is proportional to y/m*/m, where m* is the 
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effective band mass. This quantity is of the order 10~ 2 to 10" 1 [ljj, 123, indicating irrelevance (to 1/f noise) of 
quantum effects at low temperatures. Furthermore, our simulation results on transport properties |24j also suggest, 
in comparison to experiments, that quantum fluctuations do not have any visible effect with respect to the strength of 
the driving forces under consideration (see also discussion in [25I l26]p . On the other hand, due to the small parameter 
K the core action of phase slips in the bulk is large (S core oc 1/K) |22| and hence the probability of phase slips which 
are proportional to e~ core becomes small. It decreases even more under a renormalization group transformation, such 
that we can neglect phase slips in our simulations. Therefore we will use the one-dimensional classical model without 
phase slips to study the current noise in CDW systems. 



MODEL 

The charge-density p(x) of an ID CDW can be expressed as p{x) = po(l + Q~ 1 d x f(x)) + p\ cos (Qx + f{x)), where 
Q — 2kp denotes the wave vector of the undistorted wave, kp the Fermi wave vector and ip(x) a slowly varying 
phase variable. p = Q/n is the mean electron density and p\ is proportional to the amplitude of the complex order 
parameter ■ The Hamiltonian of the phase field is then given by 



/ dx | \ (JL 1 ^) ~ ^ ViS ^ x ~ X ^ x Pi cos (® x + 



n= I dx i o ( —V I ~ ^ v i S ( x ~ x i> x Pi cos (Q x + ^'( a; )) + Exd x ip{x) \ , 

where c = 4^ is the clastic constant with the Fermi- velocity vf, and Vi and Xi denote the strength and the position 
of the impurity potential acting on the CDW, respectively; E is the external electric field or driving force. 

Our numerical studies are done in the weak pinning limit, i.e. when the Fukuyama-Lee length [271 ] L c = (c/V) 2 ^ 3 
is large compared to the mean impurity distance li mp . Therefore we will restrict ourselves in the following to the case 
L c 3> hmp 3> where the full Hamiltonian can be reduced to a random field XY-model: 



H = dx \ o ( IT 1 ?) -V cos (<p - a(x)) - E<p(x)} . (2) 




Here a(x) is a random phase with zero average and e 1 ' ^ a ^ x v = S(x-x'), where the overbar denotes the averaging 

~T^T) = v 2 Sij 



over disorder realizations. V is defined by (Vipi)(VjPi) = V 2 Sij (Vi = 0). The equation of motion of the (overdamped) 



CDW is given by a Langevin equation 



where 7 is a kinetic coefficient and r](x, t) a Gaussian thermal noise characterized by (n) = and (f}(x, t) r)(x' , t')) = 
2T-f5(x-x')S(t-t'). 

1/3 

The length scale L c sets an energy scale T* = (c V 2 ) = c L~^. We will rescale time by L c /jT*, temperature 
by T*, and the external field E by E* , where E* = T* / L c is of the order of the T = depinning threshold field E c . 
In the following, E denotes the rescaled and dimensionless quantity. 

Solving the discretized version of Eq. 0one can find he time dependent current j c d w {t) which is defined as [2^] 



3cdw(t) = ( — ^ — ) , (4) 



where (...)_ denotes the average over positions. 



SIMULATION 



The effect of disorder on the dynamical behavior of the one-dimensional charge density wave model J5J) at low 
temperatures was studied l2Sl with the help of the discretized version of equation © and it was found that, contrary 
to high dimensional systems |22j , the dependency of the creep velocity on the electric field is described by an analytic 
function. The current noise spectrum was not, however, explored. Following Ref. |28j . the equation of motion is 
integrated by a modified Runge-Kutta algorithm suitable for stochastic systems with periodic boundary conditions. 
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FIG. 1: Typical time dependence of the CDW current for E = E c = 0.22 (left panel) and E = 0.3 (right panel) at T = 0.1 
for one disorder realization. The time and disorder averaged values of j c dw(t) are shown next to the curves ({jcdw))- We 
took N = 5000 and the results are averaged over 1000 and 500 samples for E = E c and E = 0.3, respectively. An initial 
(dimcnsionless) time interval of length » 2000 is discarded, such that the system is in the steady state at time 0. 



Throughout this paper, we use a system size of N = 5000 and average the results over typically N a = 1000 disorder 
realizations. Larger system sizes do not change the results substantially. Fig. 2] shows the typical time evolution of 
jcdw{t) for E = E c = 0.22 [2^| (upper panel) and E = 0.3 (lower panel) at temperature T = 0.1. One can see that the 
current exhibits strong fluctuations. The time averaged values are (j c dw) = 0.008±0.003 and (jcdw) = 0.112±0.006 for 
E = E r and E = 0.3, respectively. The spike structure is also seen, but less pronounced compared to the experimental 
data |l2j. Nevertheless, the patterns for the two values of E look similar. 

Zaitscv-Zotov studied ll2| the current noise spectrum for applied electric fields with averaged driving current 
(I) > 220pA. Using Fig. 1 from Rcf. 0] one can see that the threshold electric field in these experiments is 
E c w 35V/cm and the averaged currents of (J) = 220pA and (I) = 2.4nA at T = 2.4K correspond to electric fields 
E w 40V/cm and E w 50V/cm, respectively, i.e. the electric fields used are greater than the threshold field. Therefore 
we will restrict our spectrum analysis to E > E c . 



SCALING EXPONENT ESTIMATION 



It should be noted that in the case of non-stationary behavior of the CDW current as in our simulations, the 
standard Fourier transformation is not suitable for determining the exponent 7 and one should, therefore, employ 
more sophisticated methods. We have chosen the WTMM method |l ol| for its superior properties in non-parametric 
scaling exponent estimation |3fjj| in the presence of polynomial non-stationarities. In particular, attempts to reduce 
the non-stationary behavior of the current by discarding an initial time interval (as in^) cannot generally guarantee 
reaching a steady state, since the relaxation time to a steady state can be very long (see the remark in Ref. Q). 

The ability of the wavelet transform to provide unbiased scaling estimates of non-stationary signals is due to the 
property of orthogonality to polynomials up to the degree n of the base functions, of the so-called analyzing wavelets 
ip with to 'vanishing moments': 




x n ip(x) dx = Vn, < n < m . 



The transform is defined as the inner product of the function f(x) and the dilated and translated wavelet ip(x): 



(Wf)(s,b) = - [ dx f{x) ) , (5) 

s J s 

where s, b £ K and s > for the continuous version (CWT), which among other properties ensures local blindness to 
the polynomial bias. Indeed, the wavelet transform decomposes the signal into scale (and thus frequency) dependent 
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FIG. 2: Left: M(s) versus Ins for the CDW current j c dw{t) averaged over iV s = 1000 disorder realizations and for three values 
of the external electric field: E = 0.25, 0.3 and 0.35 which are higher than E c m 0.22. The dotted straight lines denote the 
reference slope corresponding to 7 = 1.2. The dashed line has the slope —0.5, which corresponds to the flat power spectrum 
of white noise 7 = 0. Right: The same but for T — 0.3. The dotted lines denote the reference slope H = corresponding to 
7 = 1.0. Again, the dashed line corresponds to white noise. 



components (scale and position localized wavelets) , comparable to frequency localized sines and cosines based Fourier 
decomposition, but with added position localization. This localization in both space and frequency, together with 
the wavelet's orthogonality to polynomial bias, makes it possible to access even weak scaling behavior of singularities 
h{xa)i otherwise masked by the stronger polynomial components: 

/(»)*„ =c + c 1 (x-x ) + --- + c n (x - x ) m + C\x - x \ h ^ , 

where function / is represented through its Taylor expansion around x — xq. 

In the generic multifractal formulation of the WTMM formalism jl^ , the moments q of the measure distributed on 
the WTMM tree are taken to obtain the dependency of the scaling function r(q) on the moments q: 



where Z(s,q) is the partition function of the g-th moment of the measure distributed over the wavelet transform 
maxima at the scale s considered: 



Z(s,q) = Y,(WM(s)y , (6) 

Q(s) 

with Q(s) = {uji(s)} as the set of maxima w,(s) at the scale s of the continuous wavelet transform Wf(s,t) of the 
function f(t), in our case the CDW current; f(t) = j c dwit). 

In particular, scaling analysis with WTMM is capable of revealing the modal exponent h(q = 0) for which the 
spectrum reaches maximum value; this h(q = 0) corresponds to the Hurst exponent H in the case of monofractal 
noise. This exponent is directly linked to the power spectrum exponent of the (stationary) fluctuations of the analyzed 
signal by: 7 = 2H + 1, the relation which links the spectral exponent 7 with the Hurst exponent H. 

In Fig. |2 the modal scaling exponent has been obtained by a linear fit over an appropriate scaling range from a 
suitably defined, weighted measure M{s) on the WTMM: 



h(q = 0) 



dr(q) 



dq 



q=0 



hm - — -- 

s^o log(s) 



(7) 
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FIG. 3: The dependence of the exponent h(q = 0) on temperature. Note the convergence towards H — 0.1 corresponding with 
7 = 1.2 spectral exponent for values of temperature and current. 



with 

.2(8,0) 

and for three electric field values E=0.25, 0.3 and 0.35, and for the number of the disorder averaging ensemble fixed to 
N s = 1000. Consistent with the experimental findings |l2j . the flicker noise region becomes narrower with decreasing 
E. More importantly, we obtain 7 w 1.2 as observed in experiments [l^. Asymptotic transition to the scaling regime 
characteristic to uncorrelated behavior (white noise, i.e., 7 = 0) can be clearly identified for all the values of E shown 
(see dashed line in Fig. EJ. 

Fig. [21 (right) shows M(s) versus Ins for j c dw{t) for three values of the external electric field E = 0.25, 0.3 and 0.35 
and T = 0.3. Our fitting gives 7=1, which is important from the point of view of the exact definition of l//-noisc. 
In Fig. 13 we provide the dependence of the exponent 7 on temperature. Note the convergence towards 7 = 1.2 as 
the temperature (and the averaged current) approaches 0. The exponent 7 decays quickly with temperature and we 
have the uncorrelated noise value 7 = at high T. 

The primary question remaining is that of the origins of 1// noise in the CDW system. In our opinion, the disorder 
causes the rugged energy landscape (similar to the spin glass case) leading to a wide spectrum of relaxation times. 
The average over such a spectrum would give rise to the flicker noise 0|. Our results shown in Fig. [3] support this 
point of view. Namely, at low temperatures the roughness of the energy landscape becomes more important and 
consequently the flicker-like regime occurs. Another qualitative scenario |(j for the appearance of the flicker noise in 
our system is that the CDW may be viewed as a single particle in a quasi-periodic potential with troughs of variable 
depths. Such a simplified model closely resembles the "many-pendula" model of the self-organized criticality [^H in 
which the 1/ / noise should occur. 



CONCLUSION 



In conclusion, using the classical one-dimensional CDW model and multifractal analysis, we have reproduced the 
experimental results on the current noise spectrum. Our simulations support the existence of l//-noise in this system. 
To the best of our knowledge, this is the first evidence of 1/ f scaling obtained from first principle based simulation 
in a physical (i.e. CDW) system. 

It would be interesting to check if a three-dimensional version of our model gives 7 ~ 0.8 obtained for the bulk 
NbSe3 sample ^lj or if other models should be implemented to reproduce this experimental result. The effect of higher 



6 



dimensionalities, phase slips and quantum fluctuations on the flicker noise remains a challenge for future studies. 
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